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ABSTRACT 

We exploit long-baseline ALMA sub-mm observations of the lensed star-forming galaxy SDP 81 at 
z = 3.042 to investigate the properties of inter-stellar medium on scales of 50-100 pc. The kinematics of 
the ®^CO gas within this system are well described by a rotationally-supported disk with an inclination- 
corrected rotation speed, v^t = 320 ± 20 kms“® and a dynamical mass of Mdyn = (3.5 ± 1.0) x 10®° M© 
within a radius of 1.5 kpc. The disk is gas rich and unstable, with a Toomre parameter, Q = 0.30 ± 0.10 
and so collapse in to star-forming regions with Jeans length Lj ~ 130 pc. We identify five star-forming 
regions within the ISM on these scales and show that their scaling relations between luminosity, 
line-widths and sizes are significantly offset from those typical of molecular clouds in local Galaxies 
(Larson’s relations). These offsets are likely to be caused by the high external hydrostatic pressure 
for the interstellar medium (ISM), Ptot / ^ lO^’Kcm”^, which is ^ 10"® x higher than the 

typical ISM pressure in the Milky Way. The physical conditions of the star-forming ISM and giant 
molecular clouds appears to be similar to the those found in the densest environments in the local 
Universe, such as those in the Galactic center. 

Subject headings: galaxies: starburst; galaxies: evolution; galaxies: high-redshift 


1. INTRODUCTION 

Giant molecular clouds in local galaxies follow well 
known scaling relations between CO velocity line width 
(tr) and their physical extent, P, with a oc P®/^ and mean 
moleculargas density scales and size; (n(H 2 )) oc P“® 
(|Larsonl II98II : iBolatto et al.11200811 . These scalings re¬ 
flect the dynamical state of the turbulent molecular gas 
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in the inter-stellar medium (ISM). Since most of the 
stars in local, massive spheroids and elliptical galaxies 
appear to have formed early in the history of the Uni¬ 
verse (z ~2-3), examining the physical, dynamical and 
thermal state of the molecular gas within the ISM of 
galaxies at this epoch acquires special importance. How¬ 
ever, to measure the properties of individual star-forming 
regions requires a spatial resolution of at least ~ 100 pc 
(sufficient to resolve sizes and velocity dispersions of the 
most massive giant molecular clouds). To date, this 
has only been achieved in a few rare examples of high- 
redshift galaxies whose images have be en gravitationally 
lensed by massive galaxy clu sters (e.g. iJones et'ani2010l : 
iLivermore et al.ir2012L 1201^ . 

To test whether the scaling relations that govern the 
structure of local GMGs are valid in the dense and 
rapidly evolving ISM of high-redshift, gas-rich galaxies, 
in this Letter we exploit ALMA observations of SDP 81 
- a star-forming galaxy at z = 3.042 whose image has 
been gravitationally lensed by a factor 15.8 ±0.7x by 
a ma ssive, foreground (z ^0.299) galaxy (iDve et al.l 
I2015f) . The combination of long-baselines, together 
with gravitational lensing means we are able to re- 
solve the larges t giant molecular clouds with in the ISM 
(iSolomon et al.l 119871: iScoville fc GoodI I1989H on scales 
approaching 50 pc. We use the spatially resolved 1.0 mm 
(rest-frame 250 /im) continuum imaging to identify the 
brightest star-forming regions, and measure their sizes, 
luminosities and velocity dispersions. We use a ACDM 
cosmology w ith Hn = 72kms~® Mpc~®, 11^ = 0.27 and 
Ha = 1—Hm (|Spergel et al.ll200'^ and a Chabrier IMF. 

2. OBSERVATIONS AND SOURCE PLANE 
RECONSTRUCTION 

SDP 81 was identified from the H-AT LAS survey as 
a brig ht sub-mm source at z = 3.042 bv iNeerello et all 
(|2010f) . Optical imaging and spectroscopy also revealed 

























Star Forming Regions in SDP 81 at 2 = 3.04 from ALMA 
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Fig. 1.— Source plane images of SDP81. Top Row: Left: Source-plane 1.0mm (rest-frame 250/im image) continuum image of SDP81. 
The star-forming clumps, A—E, are identified by the solid contours. The black dashed line defines the major morphological axis and the 
white contours indicate the amplification map. The luminosity weighted amplification at 1.0 mm is 16.0 ±0.7. The white ellispes show the 
source-plane PSF at the location of the clumps (offset to the right hand size of the image for clarity). Center Left: Source-plane 2.1mm 
image with the clumps identified at 1.0 mm also highlighted. The solid line denotes the caustic curve from the best-fit lens model. Center 
Right: ^^CO(5—4) velocity field. The gas disk has an observed peak-to-peak velocity gradient of 210±10kms“^. The solid white line 
shows the major kinematic axis and the dotted line denotes the major 1.0mm morphological axis. Top Right: Beam-corrected ^^CO(5-4) 
line-of-sight velocity dispersion. Bottom Row: Left: One-dimensional profile of the 1.0 mm continuum emission extracted across the major 
morphological axis. We label the five bright star-forming regions (A-E). The best-fit profiles are shown as dashed lines. Center Left: 
One-dimensional profile of the 2.1mm intensity extracted along the major axis of the 1.0mm continuum with the ^^00(5—4) intensity 
profile overlaid for comparison (blue). Center Right: One-dimensional velocity profile of the gas disk with best-fit dynamical model overlaid. 
Right: One-dimensional velocity dispersion profile (red) and intrinsic, beam-corrected (AV / Ar) velocity dispersion profile (blue). In all of 
the lower panels we show the positions of the star-forming regions. The source-plane sub-mm morphology appears complex / clumpy, but 
with the star-forming regions embedded in a dense, rotating disk. 

nosity weighted amplification is )Ltsubmm = 15.8 ±0.7 and 
Mopt = 10.2 ± 0.5 for the sub-mm and rest-frame opti¬ 
cal emission respectively. Thus, correcting for lensing 
amplification, the observed 850 /im flux density of the 
galaxy is Ssso/jm "^1.2mJy, which is representative of 
the SMG population which has rec ently been studied 
in detail, in particular with ALMA (iHodee et ahl 120131 : 
iKarim et al.l 120131: ISimpson et al.l I2014lb However, the 
amplifications means that the average source-plane reso¬ 
lution is ^ 50-100 pc - a factor ^ 30 x hi gher than that 
so-far achieved in the non- lensed case (e.g ISimpson et al.l 
120151:H^rashi et al]l2014ll . 


the presence of a massive foreground gala xy at z = 0.299 
which lenses the background galaxy (jNegrello et al.l 
1201411 . Observations of SDP 81 with ALMA in its long 
baseline configuration (up to 15 km) were taken in 2014 
Oc tober. These ob s ervati ons and reduction are described 
by iVlahakis et aP (1201511 . Briefly, the ALMA Band? 
(1.0 mm) continuum observations have a resolution of 
U X 23mas and reach a l-cr depth of 11 /xJy/beam. At 
the redshift of the galaxy, these observations sample the 
rest-frame 250/im emission. Lower-frequency observa¬ 
tions of the ^^00(5-4) and continuum emission at 2.1 mm 
w ere also made, re aching a resolution of 56 x 50 mas. 

iDve et aH (|2015[1 construct a detailed lens model for 
the system using both the ALMA sub-m m and Hubble 
Space Telescope; (HST) imaging (see also IRvbak et al.1 
1201511 . The best-fit lens model suggests that the back¬ 
ground source comprises three dominant components, a 
dense gas disk (that lies inside the caustic and so gives 
rise to the bright sub-mm emission in the image plane) 
and two galaxy nuclei (which lie on - or just outside 
- the caustic) which are visible in the HST JiL-band 
imaging. Dye et al. interpret the complex morphology 
as a merging system in which the gas disk is a result 
of an early stage interaction. In their model, the lumi¬ 


3. ANALYSIS & DISCUSSION 

Using the lens modeling from iDve et aH (1201511 . we 
reconstruct the source-plane morphology of SDP 81 at 
1.0 mm and 2.1mm and show these in Fig. [TJ The 
dust continuum morphology appears clumpy, and we iso¬ 
late five of the brightest star-forming regions from the 
highest-resolution (1.0 mm) image (using the criteria that 
they are 5-cr above the local background), and label these 
A-E. Independently, we al so use the CLUMPFIND algo¬ 
rithm (|Williams et al.l[i994ll . which isolates clumps A-D, 
although misses clump E. However, none of our conclu- 
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sions significantly change if we include / exclude region E 
from the analysis below. 

To measure the source-plane point spread function 
(PSF) we reconstruct the beam in a grid of positions 
in the image plane and measure the PSF at each re¬ 
constructed position in the source-plane. On average, 
the source-plane PSF has a FWHM^hOpc (Fig. [T]). In 
Fig. [T] we also show the one-dimensional 1.0 mm, 2.1mm 
and ^^CO(5-4) emission profiles extracted from the ma¬ 
jor morphological axis of the source. The dust contin¬ 
uum and emission profiles are not perfectly aligned. 
However, at this resolution, offsets between the CO- 
emitting gas and continuum may be expected in regions 
with high star-formation density where the ^^CO(5-4) 
(which traces the warm and dense gas) is shock-heated, 
which efficiently raises the gas temperature and density 
through mechanical heating but does not heat the dust. 
Indeed, if the disk fragments into a number of large star 
forming regions (as we discuss below), we expect large- 
scale shocks where the complexes interact. 

In Fig. [T] we also sh ow the source-pl ane ^^CO(5-4) 
velocity field (see also iDve et all l2015fl which resem¬ 
bles a rotating system with a peak-to-peak velocity of 
210±10kms“^ within 1.5kpc. The best-fit disk model 
suggests an inclination of 0 = 40 ± 5° and so a dynami¬ 
cal mass of Mdyn = (3.5 ±1.0) x 10^° Mq within a radius 
of 1.5 kpc. The total gas mass for the disk, estimated 
from either the (amplification corrected) ^^CO(l-O) lu¬ 
minosity, or using the far-infrared SED and an appro¬ 
priat e dust-to-gas ratio is 2.7-3.9 x 10^° Mq (jPve et al.l 
1201511 . Together these indicate that the central regions 
of the disk are baryon dominated with a gas fraction of 
/gas -Mgas/Mdyn-70-90%. 

The Toomre parameter, Q characterises the stabil¬ 
ity of a disk against local axisymmetric perturbations. 
Gas rich disks with Q < 1 should fragmement and col¬ 
lapse in to star forming regions. Numerical simulations 
have also suggested that tidal interactions and dynam¬ 
ical friction should force the resulting star-forming re¬ 
gions towards the center of the galaxy, where they should 
coallesce to form a bulge (which in turn stabilises the 
disk against furth er collapse; e.g. ICeverino et al.l 120101 : 
iGenzel et al.ll200^ . The Toomre parameter is calculated 
by <5 = j TrGEgas where k = a Vmax / 7? is the epicycic 
frequency (with a = \/3), is the line of sight velocity 
dispersion: and Egas is the mass surface density of the gas 
(jToomrel 1 196411 . For SDP81, we derive Q = 0.30 ±0.10, 
which is lower than typical for the gas in local ULIRGs 
(Q —1; iDownes fc Solomoiil 1 1 9981) . and slightly lower 
than the average Toomre Q of o ther gas-rich star-fo rming 
z -2 (< Q > = 0.85±0.13; IGenzel et ai.l [2(ri4l - al¬ 
though their sample also extend down to comparable val¬ 
ues, Q = 0.18 ±0.02). 

This global instability will cause large, dense gas clouds 
to form in the molecular gas. The fragmentation of the 
gas occurs on scales of the Jeans length, Lj which can be 
estimated using the gas surface density and average gas 
velocity dispersion according to Lj =7rcr^/8GSgas. For 
SDP81, with Egas = 4± 1 X 10® Mq kpc^ and crr = 30- 
35kms“^ we estimate Lj = 130/gg® pc. If the gas frag¬ 
ments on these scales, the corresponding mass of the 
clouds that form should be cr^/G^Egas — 1-2 x 10® Mq. 

Given our source-plane resolution is at least compara¬ 


ble to the Jeans length, we investigate the properties of 
the star-forming regions on these scales. As Fig.[T]shows, 
the disk contains several bright star-forming regions (A- 
E) and we measure their size, luminosity and velocity 
dispersion (using the beam-corre cted velocity disper sion 
map for the velocity dispersions: iDavies et al.ll2011h . In 
all five regions the ^®CO(5-4) has an intrinsic velocity 
dispersions of 10-35 kms“^ (Table 1). To estimate sizes 
for the molecular gas associated with these regions, we 
use two approaches. First, we fit the sub-mm light profile 
with Gaussian profiles (centered on each star-forming re¬ 
gion; Fig. P. Second, we measure the area subtended by 
a contour which is 5-cr above the local background. The 
resulting sizes (deconvolved for the source-plane PSF) 
are FWHM = 170-310 pc, which is comparable to the ini¬ 
tial fragmentation scale. We also estimate the molecular 
gas mass in the vicinity of the star-forming regions using 
the CO (5-4) line luminosity and accounting for both 
the local background and the ^^C O(l-O) / ^^CO(5-4 ) lu¬ 
minosity ratio in velocity bins Isee lDve et al.l[2015l) . ob¬ 
taining gas masses of Mgas,dumps = 1-5 x 10® Mq. These 
masses are consistent with those derived from those us¬ 
ing their line width and sizes with Mdyn = Gcr® R/G with 
G = 5 - as approriate for a uniform density sphere - with 
Mdyn /Mco = 1.5 ±0.5). 

In Fig. we compare the velocity dispersion, luminos¬ 
ity and sizes of the star-forming regions in SDP 81 with 
similar measurements for GMCs in the local Universe 
(compiled from observati ons of the Milky-W ay disk and 
other quiescent galaxies; iBolatto et al.ll200^ . as well as 
to the properties of GMCs and star-forming regions in 
more extreme environments, including the galactic cen¬ 
ter and the ISM of SMM J2135—0102 - a star-forming 
galaxy at z = 2.32 where simil ar measurements have been 
made (|Swinbank et al.l 120111) . At a fixed size, the ve¬ 
locity dispersions of the star-forming regions in SDP 81 
appear to be ^4x larger than those of GMCs in quies¬ 
cent environments, but more similar to those found in 
the Galactic Center and other gas-rich environments. In 
the local Universe, these offsets are usually attributed 
to high external p ressures on the cloud surfaces due to 
high gas de nsities (jKeto et al.ll2nn5l : iBlitz fc Rosolows^ 
l20(Ml200^ . 

To interpret the offsets in the properties of the star¬ 
forming regions in SDP 81, we use the dynamics and sur¬ 
face density of the molecular gas to estimate the ISM 
pressure (the kinetic pressure resulting from non-ordered 
mass motions). In a rotating gas disk, the mid-plane hy¬ 
drostatic pressure is given by: 


Ptot 



gas 



( 1 ) 


where Egas, E* and Cgas, cr* are the surface density 
and velocity dispersion of the gas and stars respec- 
tively. In t he Milky Way, Ptot/^B ^ 1.4 x 10^cm“®K 
(|Elinegreenl I1989D . The stellar mass of SDP 81 es¬ 
timated from the rest-frame UV / optical photometry 
is M* = (6.6l®;g) X 10^° Mq (|Negrello et al.l 1^01411 . al¬ 
though the stars appear to be offset by ^ 1.5 kpc from 
the dust and gas emission (jPve et al.l 120151 1. and so we 
adopt M* < 3 X 10^® Mq as an upper limit on the stel¬ 
lar contribution within the gas disk. We also assume 
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Fig. 2.— The relations between size, luminosity and velocity dispersions for the star-forming regions in SDP 81 compared to those in 
the both quiescent- and more extreme- evironments of the Milky Way and other local g a laxies Left: The vel ocity dispersion versus cloud 
radius. The solid line shows the line-width-size rela tion for local GM Cs (ILarsonl 119811 : IBolatto et al.ll2008ll . The normalisation of this 
scaling relation is set by the gas pressure in the disk (IElmegreenlll989r) . The dashed line shows the local relation but offset by a factor of 
10 X. GMCs in high pressure environments, such as the Galactic Center (blue points) or in gas rich galaxies such as M 64 (grey points) 
are offsets in the sense that they have higher velocity dispersions at fixed size. Indeed, the star-forming regions in SDP 81 appear to have 
velocity dispersions that are ~ 4x higher than expected. However, it is also interesting to note that their velocity dispersions are ~2 X 
lower than the star-forming regions in SMM J2135—0102. Right The relation between gas mass and size; M(H2) — R?■ The star-forming 
regions in SDP 81 are again offset from those of local GMCs, and instead are consistent with having much higer luminosities at fixed size 
- a consequence of the high turbulent pressure of the ISM (see § 3). 


that the velocity dispersion of the gas and stars are 
comparable, Ugas / o’* 1 (although allo'w this to vary 

from 0.5-2 in the calculation below). Given the high 
gas surface density, the mid-plane hydrostatic pressure 
is high, Ptot / fee x 10^ K cm“^ and although we 

caution this value has considerable uncertainty, this pres¬ 
sure is ^ 10'^ X higher than the typical pressure in the 
Milky Way disk (lO'^Kcm”^). This pressure is also 
^ 10 X higher than inferred for the ISM in more ex¬ 
treme environ ments, such as in the Galactic Center or 
the Atennae (iRosolowskv fc Blitz! 120051 : iKeto fc MversI 
119861: IWilson et al.ll200^ However, the compact disks 
of some ULIRG s, with high gas surface densities (> 
5 X 10^ Mq DC~^. lDownes fc SolomonlflQQSD may also re¬ 
sult comparably high pressures. Finally, we note that 
the implied pressure in SDP 81 is compatible with re¬ 
cent hydro-dynamic models which suggest that the typ¬ 
ical pressure in the ISM of star-forming galaxies in the 
should increases from ^ lO^Kcm”^ at z = 0.1 to ^ 10®- 
10^ K cm~^ at z = 2, rea ching ~ 10® K cm“® in some sys¬ 
tems (iCrain et al.ll2015D . 

To relate the pressure to the properties of the star¬ 
forming regions, we use “Larson ’s relations” for tu rbu- 
lent molecular clouds. Following lElmegreenl (|1989l) . the 
velocity-dispersion-size and mass-size relations can be 
cast as: 


Cr = CTo 


.fext / 

lO^Kcm”^ 



( 2 ) 


with CTo = 1.2 km s ^ (|Larso nl ll981ll . Using Equation [2l 
a pressure of 10"*^ x that of the Milky-Way therefore sug¬ 
gests that the velocity dispersions of the clouds within 
the ISM of SDP 81 should be ^lOx those in the Milky 
Way respectively (at a fixed size). The average velocity 
dispersion of the clumps in SDP 81 is four times larger 
at a fixed size than predicted from GMCs in the Milky- 
Way, and although this is lower than the factor of lOx 
predicted for pressure induced offsets alone, we reiter¬ 
ate that this is a simplified model that nevertheless re¬ 
lates the offsets in the scaling relations for GMCs with 
the gas densities and pressures in the ISM. Finally, in 
Fig. [2] we plot the mass-radius relation for a similar 
range of quiescent- and extreme- environments in both 
the Milky Way and other nearby galaxies and plot the 
positions of star-forming regions in SDP 81, which are 
again offset to higher masees by a factor ~ 10 x (at a 
hxed size) with respect to those of GMCs in quiescent 
environments. This can also be attributed to the high 
turbulent pressure (Equation [S]) . However, it is interest¬ 
ing to note that they do not appear as massive as the 
star-forming regions in SMM J2135—0102, although the 
latter is a natural consequence of the mass scale for col¬ 
lapse within the disk given their relative velocity disper¬ 
sions and gas surface densities (ctsmm J 2135 / crsopsi 2.3 
and Egas,SMMj2i35/Sgas,sDP8i ^2) which results 
in a mass scale difference for the clumps of 
Mci,smmj 2135 / Mci,sdp 81 ^ 25 - Consistent with Fig. [TJ 

4. GONCLUSIONS 


and 


Pext/fcB 
10^ cm“3 K / 



(3) 


Using long-baseline ALMA observations, we have 
mapped the distribution of star-formation and molec¬ 
ular gas in the lensed, star-forming z = 3.042 galaxy, 
SDP 81, on physical scales of ~ 50-100pc. The ^®CO(5- 


M(H2) = 290 
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TABLE 1 
Clump Properties 


ID 

^clump 

(km 

Amplification 

(m) 

FWHMalamp 

(pc) 

^clump 

(km s~^) 

/cOC5-4) 
(mjy kms“^) 

r54 

Mgas 

(x 10^ Mo) 

A 

-42 ± 65 

7.2 ±0.8 

282 ± 25 

30 ±9 

57 ±4 

0.28 ±0.05 

1.2 

B 

162 ±12 

40.1 ±1.0 

188 ± 25 

35 ±3 

24 ±4 

0.30 ±0.12 

2.8 

C 

-150 ±28 

6.2±0.18 

305 ± 95 

31±7 

46 ±4 

0.28 ±0.08 

4.8 

D 

-176 ±2 

5.9±0.15 

300 ± 70 

11±7 

35 ±4 

0.36 ±0.05 

0.8 

E 

204 ±5 

30.4 ±7.4 

170 ±40 

31 ±3 

28 ±4 

0.18 ±0.04 

0.7 


Notes: i^clump denotes the velocity of the disk at the position of the star-forming region with respect to the systemic redshift. The 
amplifications are 1.0 mm emission luminosity weighted values. FWHMciump have been deconvolved for the source-plane PSF. (rdump is 
the velocity dispersion of the clump as measured from the intrinsic ^^CO(5-4) velocity dispersion map. is the ratio of the ^^CO(5- 
4) /^^CO(l-O) luminosities derived from velocity and amplification maps; JDve et al.ll201^ . Mgas denotes the gas mass of each clump 
assuming Mgas = aco Lco(l-O) Lco(i_o) = Lco(i-o) / ''54 and “CO(l-O) = Mgas / L'co(l-o) = 


4) dynamics suggest that the molecular gas is lo¬ 
cated in a disk with an inclination-corrected rotation 
speed of lij-ot =320±20kms“^ and a dynamical mass of 
Mdyn = (3.5 ±1.0) X 10^° M© within a radius of 1.5 kpc. 

The gas disk appears to be Toomre unstable, 
Q = 0.30 ± 0.10. This instability will cause large, dense 
star-forming regions to collapse on scales of the Jeans 
length, Lj = 130;gQ° pc. We identify five star-forming 
regions on these scales in the rest-frame 250-^m contin¬ 
uum and measure their size, luminosity and CO velocity 
dispersions. We show that these star-forming regions do 
not lie on the local relations for GMCs, but are instead 
systematically offset such that the velocity dispersion is 
^ 4 X higher than typical GMCs at a fixed size. 

The gas dynamics and surface density suggest that the 
ISM should be highly pressurised, and we estimate a mid¬ 
plane hydrostatic pressure which is ^ 10^ x higher than 
typically found in the Milky Way. These high pressures 
appear to be responsible for the in offsets in the scaling 
relations of the star-forming regions compared to those of 
typical GMCs in galaxies in the local Universe. Within 


the star-forming ISM of this dense gas disk, the phys¬ 
ical conditions appear to be similar to those only seen 
in the densest environments in the local Universe (e.g. 
'Kruiissen fc Longmorell20i^ 
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